Effect of frustration on charge dynamics for a doped two-dimensional triangular 
Hubbard lattice: Comparison with a square lattice 
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We examine the optical conductivity a(ui) and the chemical potential fi, together with the spin 
correlation, in the strong-coupling limit of a hole-doped two-dimensional triangular Hubbard model 
near half filling by using an exact diagonalization technique. In contrast to the case of a square lattice 
without frustration, the doping dependences of fj, and the Drude weight indicate that the charge 
degree of freedom is weakly coupled to the spin degree of freedom. However, we find that cr(u>) shows 
strong incoherent excitations extended to a higher energy region. This implies that geometrical 
frustration in strongly correlated electron systems influences incoherent charge dynamics. 
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The interplay of geometrical frustration and strong 
electron correlation is one of the hot topics in the field 
of strongly correlated electron systems. The simplest ex- 
ample of such geometry is a triangular lattice in two di- 
mensions (2D)i 

Numerical studies for a spin- 1/2 triangular Heisenberg 
model have suggested the presence of the antiferromag- 
netic (AF) long-range order of a three-sublattice copla- 
nar typ o 2 '?' in spite of frustration. However, the order 
is expected to be easily destroyed by carrier doping and 
a resonant-valance-bond (RVB) stated may emerge. In 
fact, it has been suggested that for the positive sign of 
hopping (t > 0) the RVB state is favored with hole dop- 
ing& and d + id-wave superconductivity appears^ It was 
also suggested that the three-sublattice magnetism is sta- 
ble in a wide range of doping^ On the other hand, for 
t < the Nagaoka ferromagnetism^ emerges with dop- 
j n g | 6 J 7 J 8 i io These results indicate that geometrical frus- 
tration and strong correlation give rise to strong compe- 
tition among many states. 

The variety of the ground states in the triangular lat- 
tice may induce unconventional electron excitations and 
charge dynamics through coupling with frustrated mag- 
netism. However, there are few theoretical studies about 
charge dynamics away from half filling, except for the 
fermion-spin theory of the t- J modeliii*^ Therefore, it is 
very important to investigate the effect of frustration on 
charge dynamics and to clarify the interplay of the charge 
and spin degrees of freedom in geometrically frustrated 
systems. For this purpose, comparison with a square 
lattice is very crucial, because it has no geometrical frus- 
tration. 

In this paper, we perform an exact diagonalization 
study of the doping dependences of the optical conduc- 
tivity <j(u), the chemical potential n, and static spin cor- 
relation for a twenty-one-site triangular Hubbard clus- 
ter with large on-site Coulomb interaction. We intro- 
duce an averaging procedure over various twisted bound- 
ary conditions (BCs) to reduce finite-size effects. From 
spin correlation near half filling, we confirm strong three- 
sublattice magnetic correlation for t > and a ten- 
dency toward ferromagnetism for t < 0, as reported be- 



f ore> 6 i 7 i 8 4 io w e fj nc j -that /i and the Drude weight D are 
roughly proportional to the hole concentration 6, irre- 
spective of the sign of t. These behaviors are similar to 
those of the one-dimensional (ID) Hubbard model where 
the charge and spin degrees of freedom are decoupled, 
but different from the case of the square lattice where 
the coupling is strong. Therefore, charge in the trian- 
gular lattice is expected to be weakly coupled to spin. 
However, we find that cr(w) shows very strong incoherent 
excitations extended to a high energy region of several t. 
Such high-energy incoherent excitations are not seen in 
the square lattice. This implies that geometrical frustra- 
tion causes unconventional charge dynamics. 

The Hubbard model with the nearest-neighbor hop- 
ping t and on-site Coulomb interaction U is given by 

#Hub = -t (ctffCj'^ + H.C.J + U^2ni,irii,i, (1) 

where c^ a annihilates an electron with spin a at site 
i> n-i,a — cl a Ci^, and the summation runs over 

nearest-neighbor pairs. Being interested in the region 
of U ^> t, we take the strong-coupling limit of Eq. QJ. 
The resulting Hamiltonian reads Hsc — HtJ + H3S with 




(i,j,k),<r,a' 



Xrij-e'Cy ^Ck,*' (1 - nk,-a>) , 

where J = U 2 /U, c iyCr = Ci t< ,{l-TH,- a ), hi = ct 
and (i, j, k) denotes a pair of three nearest-neighbor sites. 
We use U/t = 20, which is above a critical U/t [~ 12.1 
(Ref. 13) or ~12.4 (Ref. 14] for the metal-insulator (MI) 
transition at half filling. 

The exact diagonalization method based on the Lanc- 
zos algorithm is frequently applied to the Hubbard-type 
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models as an unbiased numerical method. For the trian- 
gular lattice, one uses a TV-site cluster with R a = lu+mv 
and R;, = -mu + (I + m) v, being that N = I 2 + rn 2 + Im 
with integers l,m > and u and v are the vectors 
connecting nearest-neighbor sites given by u = x and 
v = -^y with the unit vector x (y) in the x (y) di- 
rection. In this study, we take TV = 21 (I — 4, m — 1) and 
introduce carriers up to five holes. In such a small cluster, 
we are not free from finite-size effects. In order to reduce 
the finite-size effects, we introduce various BC with twist 
and average physical quantities over the twisted BC. This 
procedure has been applied for various quantities in t-J- 
type modelsii^i^ The twist induces the condition that 
Ct+R a ,a = ^ a Ci >a and c i+Rb , a = e'^c^ with arbitrary 
phases 4> a and (/){,. Note that 4> a = 4>b = corresponds 
to periodic BC. 4> a (b) is defined as <p a (b) = K ' R-a(h) with 
hi/ — /s^x "T" fcyy . h\i is taken within an area surrounded 
by four corners at (k x ,k v ) — ±j^(l,^(l + 2m)) and 

±%(l + 2m,-V30- For the averaging procedure, we 
choose many n with equal intervals of 7r/40 in the area. 
The total number of k, results in TV K = 353, except for the 
five-hole case where we take TV K = 177 because of time- 
consuming calculations. For the square lattice, we use 
a TV = 20 cluster^ and TV K = 320. We note that, at 
half filling (5=0) where the Hamiltonian is given by the 
Heisenberg model, the averaging procedure can not be 
applied because of the absence of hopping terms. 

We first examine the static spin correlation function 
C s (r) between two sites separated by a distance r. Using 
the averaging procedure, we define this as 
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where I^q) represents the ground state with energy E* 
for a k, the summation of p runs over bonds satisfying 
\p\ = r, and N r is the number of the bonds. Figure 
shows the doping dependence of C s (r). In the triangu- 
lar lattice with t > [Figs. GJa) and[Hb)], C s (r) for 
5 > shows r dependences similar to that for 5 = 0. 
This behavior is consistent with a previous report^ that 
the three-sublattice magnetism at half filling is stable in 
a wide range of 5. We also find that this doping de- 
pendence is different from the case of the square lattice 
shown in Figs.[TJe) and Off), where AF spin correlation 
at 5 = is rapidly suppressed with doping because of the 
strong coupling of charge motion and spin background. 
Therefore, this comparison may tempt us to have a view 
that the spin and charge degrees of freedom are weakly 
coupled in the t > triangular lattice. However, such a 
view is different from the spin-charge separation seen in 
ID systems, as will be discussed below in terms of <j(uj). 
The t < triangular lattice shows a different behavior. 
Even for single- hole doping, C s (r) dramatically changes 
from that of the Heisenberg model as shown in Fig.^c), 
and approaches to small values toward 5 = 5/21. This 
may be a precursor of ferromagnetic ground states i 6 i 7,8 i 10 
In fact, the total spin at 5 = 5/21 under the averaging 
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FIG. 1: Spin correlation C s as a function of two-spin distance 
r for various 5 in the strong-coupling limit of the Hubbard 
model with U/t = 20. (a),(b) the N = 21 t > triangular 
lattice, (c), (d) the N = 21 1 < triangular lattice. (e),(f) the 
TV = 20 square lattice. The numbers in the panels represent 
5. 



procedure is = 0.441. We note that a fully polar- 
ized ferromagnetic ground state is obtained at 5 = 7/21 
under periodic BC. It is interesting to notice that the 
doping dependence of C s (r) is small except for the dra- 
matic change from zero to one hole. This is also different 
from the case of the square lattice, where the change from 
zero to three holes is continuous. The small doping de- 
pendence for t < also seems to be indicative of weak 
coupling of spin and charge. 

In small clusters, the chemical potential fi at a given 
5 can be expressed as \i = [eo(<5i) — eo(<$2)] l(5i — 61), 
where eo(5) is the ground-state energy per site evalu- 
ated by e (S) = (N k N)- 1 Y ik Eq and 8 X and 5 2 sat- 
isfy 5 = (Si + <5 2 ) / 2 - Here we note that the averag- 
ing for eg(5) is crucial for obtaining reliable doping de- 
pendence of fi because the averaging reduces the finite- 
size effects and gives ground-state energies that smoothly 
change with carrier concentration. By using Eq , the 
Drude weight D along the x direction is given by D = 
N- 1 J2 K (2N)-id 2 E£/d4- 

In Fig.0 [i is plotted as a function of 5. /1 for both the 
t > and t < triangular lattices follows well the linear 
solid lines, implying that fi is proportional to 5. On the 
other hand, /1 in the square lattice clearly deviates from 
a linear relation, but is well described by a quadratic 
polynomial. The relation that \i oc S 2 is consistent with 
previous numerical resultsiii^ This contrasting behav- 
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FIG. 2: Chemical potential /j, versus <5 in the strong-coupling 
limit of the Hubbard model with U/t = 20. The upper and 
lower triangles represent the data for the N = 21 triangular 
lattice with t > and t < 0, respectively, and the squares are 
for the N = 20 square lattice. The filled and open symbols 
correspond to /i obtained by taking the concentration differ- 
ence of S-2 — Si = l/N and 2/N, respectively. The solid lines 
are guides for eyes with linear and quadratic functions for the 
triangular and square lattices, respectively. 



ior between the triangular and square lattices indicates 
essential difference of their electronic states. 

According to quantum criticality hypothesis of the MI 
transition, a hyper scaling theory has been developed*^ 
The dynamical critical exponent z is related to the doping 
dependence of fi in the form that fx — [i c cx S z / d , where d 
is the space dimension and fi c is the chemical potential at 
a critical concentration (6 = 0). The relation that z = 2 
is obtained for the ID Hubbard model from the exact 
solution^ From Fig. |2 the d — 2 triangular lattice also 
satisfies z = 2. However, in the square lattice, z = 4 
since [i — [i c cx S 2 i& 

The hyperscaling theory also predicts that D cx 
( ^i+(z-2)/d ) i9 g mce z — 2 f or the triangular lattices, we 
expect that D cx 8. Figures |3Ja) and|3Jb) show the dop- 
ing dependece of D. The open and solid circles represent 
D obtained for the periodic BC and for the averaging 
procedure, respectively. Although D does not necessar- 
ily show smooth behaviors because of the finite-size ef- 
fect, it seems to grow roughly proportional to 8 with a 
slight convex behavior below 8 ~ 0.2, i.e., D cx S a with 
a ^ 1, for both the t > and t < cases. It is diffi- 
cult to estimate the value of a from these data, but we 
consider that the data do not contradict the predicted 
behavior of D cx 8 from /i. This suggests that the scal- 
ing theory is satisfied in the triangular Hubbard model. 
Since the scaling relation with z — 2 is also satisfied in 
the ID Hubbard model where charge is decoupled from 
spin in the ground state, we can say that, even in the 
triangular lattice, charge would be coupled very weakly 
to spin. This view is consistent with the results of spin 
correlation discussed above. 

In the square lattice, two different behaviors have 
been reported: D is proportional to 8 2 (Ref. 21) and 
8 (Refs. 22,23). Figure EJc) shows proportionality to 8. 
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FIG. 3: Drude weight D versus 8 in the strong-coupling limit 
of the Hubbard model with U/t = 20 for the N = 21 trian- 
gular lattice with (a) t > 0, (b) t < 0, and (c) the N = 20 
square lattice. The solid circles are obtained by averaging 
over twisted BC. The open circles represent the data under 
periodic BC and for the even number of electrons, (d) Dop- 
ing dependence of the integrated weight of a Teg (u>). The upper 
and lower triangles represent the data for the N = 21 triangu- 
lar lattice with t > and t < 0, respectively, and the squares 
are for the N = 20 square lattice. 



The resulting z = 2 is different from z — 4 deduced from 
fi. This inconsistency may tell us that the MI transition 
of the square Hubbard lattice is not a continuous phase 
transition. This is in contrast with the ID case and prob- 
ably with the triangular lattice. The difference might be 
attributed to the presence of strong coupling of charge 
and spin in the square lattice. However, we have to say 
that our numerical data is still not enough to give a fi- 
nal conclusion, because D very close to half filling is not 
rigorously linear in 8. Therefore, we need further careful 
efforts to solve this important problem. 

The real part of the optical conductivity reads cr(uj) = 
2ttD8(uj) + cr rcg ((jj), where the regular part <7 reg (w) under 
the averaging procedure is given by 



where |\&5n) represents an eigenstate with energy for 
a given k. The x component of the current operator is 
given by j x = i [Hsc j x ] i where x is the x component 
of the total position operator. A standard continued- 
fraction expansion method based on the Lanczos algo- 
rithm is used to calculate a TCg (uj). 

Figure I3J shows the doping dependence of a(oj). In the 
triangular lattices, we find that incoherent excitations are 
widely extended to a high-energy region up to several t. 
This sounds counterintuitive from the viewpoint of weak 
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FIG. 4: (Color online) Optical conductivity o{ui) for various 
hole concentrations of the strong-coupling limit of the Hub- 
bard model with U/t = 20. (a) t > 0, (b) t < for the 
TV = 21 triangular lattice, and (c) the TV = 20 square lattice. 
The data are obtained with the Lorentzian broadening of 0.1/. 
The numbers in the panels represent the hole concentration. 



coupling of charge to spin in the ground state mentioned 
above, since, in the ID Hubbard model where spin and 
charge are completely decoupled, incoherent weights are 
negligibly small^4 The incoherent excitations are also dif- 
ferent from those of the square lattice. For the square lat- 
tice, the origin of the incoherent part is well-understood 
based on the picture that magnetic polarons move in the 
system by inducing magnetic excitations scaled by the 
exchange interaction J = 4t 2 /U ~ 0.2t. Therefore, the 
incoherent spectral weight shown in Fig. |3J C ) is mainly 
located below u *~ t. Very incoherent <j{u>) up to several 
t in the triangular lattices clearly indicates that such a 
magnetic-polaron picture cannot be applied. 



An interesting feature of <j{u>) for both the t > and 
t < triangular lattices is that the incoherent spectral 
weight increases with increasing 6 without changing their 
global energy distribution in the region of t < lo < 5t. 
Furthermore, the spectral weight of the incoherent part 
is larger than that of the square lattice, as shown in 
Fig- Eld) where the integrated weight of a Ieg (u>) is plot- 
ted. Since the magnitude of D is comparable between 
the triangular and square lattices, the large incoherent 
weight in the triangular lattice means that the motion of 
carrier is more incoherent than the square lattice. En- 
hanced incoherent structures have also discussed in terms 
of the single-hole spectral function in triangular antifer- 
romagnet:— Both quantities may have common origins. 
A possible origin is nontrivial phases from the spin de- 
grees of freedom^ with strong frustration. The phases 
may act as random fields on carriers and induce inco- 
herent excitations. Obtaining analytical expressions for 
theses effects is a future problem. We note that such 
an enhanced incohcrcncy might be a possible origin of 
anomalous incoherent contribution of <r(ui) observed in a 
triangular compound Nao.7Co02»2£ 

In summary, we have performed the exact diagonal- 
ization study of the doping dependence of the optical 
conductivity cr(u>), the chemical potential /i, and static 
spin correlation for the triangular Hubbard model. We 
find that /i and the Drude weight D near half filling are 
roughly proportional to the hole concentration. These 
behaviors are similar to those of the ID Hubbard model 
where the charge degree of freedom is decoupled from 
spin, but different from the case of the square lattice 
where the coupling is expected to be strong. Therefore, 
charge in the triangular lattice is expected to be weakly 
coupled to spin. Spin correlation function near half filling 
supports this view. However, we find that ct(cj) shows 
very strong incoherent excitations, which are absent in 
the square lattice. This implies that geometrical frus- 
tration influences incoherent charge dynamics. A clear 
understanding of the physics behind this remains for the 
future. 
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